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August, 2000 
Abstract 

From a sample of 848 it 44 — > 7r^7r+7r+ decays, we find T{Df — > 
7r-7r+7r+)/r(Z)+ ^ <j)TT+) = 0.245 ± 0.028/[!;[ji^. Using a Dalitz plot analysis 
of this three body decay, we find significant contributions from the channels 
p°(770)7r+, p0(1450)7r+, /o(980)7r+, /2(1270)7r+, and /o(1370)7r+. We present 
also the values obtained for masses and widths of the resonances /o(980) and 
/o(1370). 
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The charm meson decay ti^h^h^ and its charge conjugate (imphcit 

throughout this paper) is Cabibbo-favored but has no strange meson in the final 
state. The decay can proceed via spectator amphtudes, producing intermediate 
resonant states with hidden strangeness, e.g. —>■ /Q(980)7r"'", with ss quarks in 
the /o(980). Also, the decay can proceed via VF-annihilation amplitudes, produc- 
ing intermediate states with no strangeness, e.g. —>■ p^n^. A H^-annihilation 
amplitude could also produce the intermediate state /Q(1370)7r+, assuming 

the /o(1370) consists mostly of uu and dd quarks as predicted by the simple quark 
modelJU. To determine the relative importance of these different decay mechanisms, 
one can use an amplitude analysis. Such an analysis is also able to determine the 
masses and decay widths of the intermediate states. 

In general, scalar resonances have large decay fractions in the 3-body decays 
of Z)-mesons, and such decays provide a relatively clean laboratory in which to 
study the properties of the scalars. In particular, isoscalar intermediate states are 
dominant in 7r^7r+7r+ decays[0, 0. The largest contribution to this final state 

comes from the decay involving the scalar meson /o(980) whose nature is a long- 
standing puzzle. It has been described as a qq state , a KK molecule, a glueball, 
and a 4-quark statef^. 

In this paper we extend the reach of previous studies using the larger data sample 
from Fermilab experiment E791. We present an amplitude analysis which includes a 
greater number of possible resonant states, and we measure the masses and widths 
of the scalar resonances /o(980) and /o(1370) with better precision. Taken together 
with the results of the companion analysis[^], these results provide new insights into 
the characteristics of scalar mesons and their importance in charm meson decay. 

The data were produced by 500 GeV/c vr" interactions in five thin foils (one 
platinum, four diamond) separated by gaps of 1.34 to 1.39 cm. The detector, the 
data set, the reconstruction, and the resulting vertex resolutions have been described 
previously 0. After reconstruction, events with evidence of well-separated produc- 
tion (primary) and decay (secondary) vertices were retained for further analysis. 
From the 3-prong secondary vertex candidates, we select a 7r~7r~'"7r+ sample with in- 
variant mass ranging from 1.7 to 2.1 GeV/c^. For this analysis all charged particles 
are taken to be pions; i.e., no direct use is made of particle identification. 

We require a candidate's secondary vertex position to be cleanly separated from 
the event's primary vertex position and from the closest target material. The sum 
of the momentum vectors of the three tracks from this secondary vertex must point 
to the primary vertex. The candidate's daughter tracks must pass closer to the 
secondary vertex than to the primary vertex, and must not point back to the primary 
vertex. The resulting invariant mass spectrum is shown in Figure |I|. 

We fit the spectrum of Figure |1| as the sum of and signals plus back- 
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ground. To account for the signal's non Gaussian tails, we model each signal as 
the sum of two Gaussian distributions with the same centroidbut different widths. 
We model the background as the sum of four components: a general combinatorial 
background, the reflection of the — > K~tt^t[^ decay, reflections of — > K~t:^ 
plus one extra track (mostly from the primary vertex), and — > t/'tt^ followed 
by 7]' p'^(770)7, p°(770) tt+tt^. The K^-k^-k^ reflection is located be- 

low 1.85 GeV/c^ in the 7r^7r+7r+ spectrum. The other charm backgrounds populate 
the whole tt~t[^tt^ spectrum. We use Monte Carlo (MC) simulations to determine 
the shape of each identified charm background in the tt'tt^h^ spectrum. We as- 
sume that the combinatorial background falls exponentially with mass. The levels of 

— s> K^Ti^ and rj'-K^ backgrounds are determined using charm signal rates 

measured in our total event sample and branching ratios taken from the compilation 
by the Particle Data Group The parameters describing the combinatorial back- 
ground and the level of the K^ti^h^ reflection are determined from fitting 
the 7r^7r+7r+ distribution. The mass (centroid) and both Gaussian widths for each 
signal float in our fit. The fit finds 1172 ± 61 events and 848 ± 44 events. 
The Z)+ signal is used for the measurement of the decay rate T{D^ 7r^7r+7r+) 
relative to T{D^ — > ^vr"*") and for an amplitude analysis of the Dalitz plot. 

To minimize systematic effects when calculating the ratio of efficiencies, we select 
the — >■ 071"*", — > K~K^ signal using the same track and vertex quality criteria 
as for the — 7r^7r+7r+ decay. The number of normalization events is found to 
be 1038 ± 44. 

We use MC simulations to correct the signals for geometrical acceptance and 
detector efficiency. We measure a ratio of efficiencies e{D^ — >• 7r~7r+7r+)/e(Z}+ — >■ 
07r+) = 1.64 ± 0.15 with values e{D+ 7r-7r+7r+) and e{D+ 07r+) of about 2% 
and 1%, respectively. The error systematic is dominated by uncertainties in the MC 
model of D production and relative efficiencies of our selection criteria. 

The branching fraction for 7r~7r+7r+ relative to that for (f)7i~^ is 

measured to be: 

B(Dt 7r-7r+7r+) .^m^ 

The first error is statistical. The second is systematic, and is dominated by uncer- 
tainties related to the signal and background shapes used in the fit, the background 
levels, and the sample selection criteria. This value is smaller than the other experi- 
mental resuhs: E691 § found 0.44 ±0.10 ±0.04, WA82 [0 quoted 0.33 ±0.10 ±0.04 
and E687 § presented 0.33 ± 0.058 ± 0.058. The PDG ^ presents a value of 
0.28 ± 0.06 from a constrained fit. 

The symmetrized Dalitz plot of the 937 candidates with invariant mass between 
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Figure 1: The h'ti^h^ effective mass spectrum. The dotted hue represents the 
— > i^~7r+ plus Tj'Ti^ and the dashed hne is the total backgound. Events 

in the hatched areas at the Dg mass are used for the Dalitz plot analysis in this 
Letter. The hatched area at the mass is used for the analysis in the following 
companion paper [4]. 

1.95 and 1.99 GeV/c^ is shown in Fig. 0. The integrated signal to background 
ratio is ~ 2. The narrow horizontal and vertical bands of Si2 = 771^(71^7?^) and 
Si3 = m^(7rf7r^) just below 1 GeV^/c^ correspond to the /o(980)7r+ state with 
constructive quantum-mechanical interference evident where the two bands overlap 
(the event count is four times that of the individual bands). At the upper edge of 
the diagonal, there is another concentration of events centered at S12 ^ S13 ~1.8 
GeVVc^ corresponding to the /2(1270)7r+, /o(1370)7r+, and p°(1450)7r+. 

We fit the distribution shown in Figure ^ to a signal probability distribution 
function (PDF), which is a coherent sum of amplitudes corresponding to the non- 
resonant decay plus five different resonant channels, and a background PDF of known 
shape and magnitude. The resonant channels we include in the fit are p''(770)7r+, 
/o(980)7r+, /2(1270)7r+, /o(1370)7r+, and p0(1450)7r+. We weight the signal PDF by 
the acceptance across the Dalitz plot and use the measured line shape and back- 
ground to determine the ratio of expected signal to background for each event in 
the Dalitz plot. 

We assume the non-resonant amplitude to be uniform across the Dalitz plot. 
Each resonant amplitude, except that for the /o(980), is parameterized as a product 
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of form factors, a relativistic Breit-Wigner function, and an angular momentum 
amplitude which depends on the spin of the resonance. 



A - ^ " M'' (2) 

with 

The form factors Fo and "^F^ are the Blatt-Weisskopf damping factors respec- 
tively for the D and the resonance decays, p* is the pion momentum in the resonance 
rest frame at mass m^r^ (pg = p*(mo)) and J is the spin of resonance n. M.f^ describes 
the angular distribution due to the spin[|]. Since we have identical particles in the 
final state, each signal amplitude is Bose-symmetrized, An = ^„[(12)3] +^„[(13)2]. 

For the /o(980)7r"'" we use a coup led- channel Breit-Wigner function, following the 
parameterization of the WA76 Collaboration Q, 

with 

and 



= g^Jml^/A~ml (5) 



/A-ml^ + ^mlJA-ml,)j . (6) 

We multiply each amplitude by a complex coefficient, Cj = aje^^ . The fit param- 
eters are the magnitudes, aj, and the phases, 6j, which accomodate the final state 
interactions, are fit parameters obtained using the maximum-likelihood method. 

Monte Carlo simulations are used to determine the shape and location of the 

— > K~7r~^ background in the Dalitz plot. The amount the background is 
determined using both MC simulations and data. The contribution of the 
rj'-K^ background is negligible. The background proportions are 12 ± 2% for — > 
K^TT^ and 88 ± 2% for the combinatoric across the Dalitz plot. Checks of the 
background model are described in the companion paper ||^. 

The parameters of the /o(980) state, g.,,, gx, and mo, as well as the mass and 
width of the /o(1370), are determined directly from the data, fioating them as free 
parameters in the fit. The other resonance masses and widths are taken from the 
PDG|T|. The results of the D'^ Dalitz plot fit are shown in Table I. The column 
labeled Fit A corresponds to our best fit with all six modes. The measured values 
mo = 977±3±2 MeV/c^, g^ = 0.09 ± 0.01 ± 0.01 and gx = 0.02 ± 0.04 ± 0.03 have 
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Figure 2: The Df — >■ vr vr+vr"'" Dalitz plot. Since there are two identical particles, 
the plot is symmetrized. 

been corrected for small shifts in these parameters due to the 7r^7r+ mass resolution 
in this region. This resolution is estimated to be approximately 9 MeV/c^. We 
determine the shifts by folding a 9 MeV/c^ resolution together with various qk, Qn 
combinations and fitting the resulting distributions to the form in Eqns. (4)- (6). 
The resolution was found to shift mo, Qk^ and g.„ by —1.4 MeV, +0.06 and —0.006 
respectively. These shifts have been included in the values quoted above and in 
Table I. Uncertainties in our resolution contribute to the systematic uncertainties. 
The magnitudes and phases of the resonant amplitudes are relatively insensitive to 
the value of qk- These values are not compatible with WA76 results^l, Qt, = 0.28 
± 0.04 and gx = 0.56 ± 0.18. For the /o(1370) we find mo =1434 ± 18 MeV/c^ and 
Fo = 173 ± 32 MeV/c^. We have also fit the Dalitz plot using for the /o(980) the 
same Breit-Wigner function as for the other resonances. The resulting fit is nearly as 
good as the one using the coupled-channel Breit-Wigner function, and the fractions 
and phases are indistinguishable. With this parameterization we find mo =975 ± 3 
MeV/c^ and Fq = 44 ± 2 ± 2 MeV/c^, also shifted by the effect of mass resolution. 

Table I shows the magnitudes {aj) and phases determined from the fit, and 
corresponding fraction for each decay mode. The fractions are defined as 



Jdsudsis 


CjAj 


2 


J dSi2dSi3 J2jk 


CjAjC 


* A* 1 
k-^k 1 



(7) 
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The first reported error is statistical and the second is systematic, the latter being 
dominated by the uncertainties in the resonance parameters, in the background 
parameterization, and in the acceptance correction. The /o(980)7r"^ is the dominant 
component, accounting for nearly half of the —>■ 7r~7r"'"7r"'" decay width, followed 
by the /o(1370)7r+ and /2(1270)7r+ components. The contribution of p°(770)7r"'" 
and p''(1450)7r+ components corresponds to about 10% of the 7r^7r+7r+ width. We 
have not found a statistically significant non-resonant component. The S12 and S13 
projections are nearly independent and the sum of the two is shown in Fig. ^ for 
Fit A. 

To assess the quality of our fit absolutely, and to compare it with other possible 
fits, we developed a fast-MC algorithm that simulates the ti^ti^h^ Dalitz 

plot from a given signal distribution, background, detector resolution and accep- 
tance. For any given set of input parameters we calculated a using the procedure 
presented in Ref. [Q. From and the number of degrees of freedom (z/), we cal- 
culate a confidence level assuming a Gaussian distribution in x^/^- The confidence 
level for the agreement of the projection of Fit A onto the Dalitz plot with the data 
of Fit A is 35%. 

We perform fits excluding amplitudes with small contributions. The fit without 
the non-resonant amplitude is as good as Fit A, and the resulting parameters are 
essentially the same. When comparing Fit A with models without the p^vr"*" ampli- 
tudes (Table 1) we calculate Aw = — 2(ln£j — IiiCa), where Ca is the likelihood of 
the fit with all modes, and is the likelihood of the different models, as described in 
the companion paper]^. For the model where we exclude the p^{770)n~^ amplitude 
we have, for data. Aw =20. In the fast-MC with all modes we have (Aw) = 37. 
In the fast-MC with no p'^(770)7r"'" we have (Aw) = —29. We observe a similar 
behavior in Fit C, where we exclude the p°(1450)7r"'~ amplitude, we have Aw = 32. 
In the fast-MC with all modes we have (Aw) = 30 while with the fast-MC with no 
p''(1450)7r+ (Aw) = —30. In all fast-MC experiments the rms deviation for (Aw) 
is about 11 units. We conclude that the best description of our data includes both 
p°(770)7r+ and p°(1450)7r+ amplitudes. 

The contribution of modes having isoscalar mesons completely dominates the 
TT^TT+TT^ decay. The same isoscalar dominance is observed in the — > 
7r^7r+7r+ decay Q. However, there is no evidence in the D'^ decay for a low-mass 
broad scalar particle as seen in the decay. If the 7r~7r"'"7r"'" decay is 

dominated by the Cabibbo-favored spectator mechanism, we would expect final 
states with a large ss content. Approximately half of the 7r^7r+7r+ rate is 

produced via /o(980)7r^. The /o(980) is often supposed to have a large ss component, 
indicating a large spectator amplitude in this decay. On the other hand, the large 
contribution from the intermediate state /o(1370)7r"'" indicates the presence of either 
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Fit A 


Fit B 


Fit C 




Fraction(%) 


Fraction(%) 


Fraction(%) 




Magnitude 


Magnitude 


Magnitude 




Phase 


Phase 


Phase 


/o(980)7r+ 


56.5 ± 4.3 ± 4.7 


58.0 ± 4.9 


54.1 ± 4.0 




1 (fixed) 


1 (fixed) 


1 (fixed) 




U(nxedj 


0(nxedj 


0(nxedj 


NR 


0.5 ± 1.4 ± 1.7 


7.5 ± 4.8 


5.0 ± 3.8 




0.09 ± 0.14 ± 0.04 


0.36 ± 0.12 


0.30 ± 0.12 




(181 ± 94 ± 51) 


(165 ± 23) 


(149 ± 25) 


pO(770)7r+ 


5.8 ± 2.3 ± 3.7 





11.1 ± 2.5 




0.32 ± 0.07 ± 0.19 





0.45 ± 0.06 




/inn _L o /I _L 
(iOy ± 24 ± 5j 





(81 ± 15j 


/2(1270)7r+ 


19.7 ± 3.3 ± 0.6 


22.2 ± 3.3 


20.8 ± 3.0 




0.59 ± 0.06 ± 0.02 


0.62 ± 0.06 


0.62 ± 0.05 




[loo ± io ± zoj 


/I no 1 1 1 \ o 

(^iuy ± iij 


/I 0/1 1 1 1 \ 
[iZ4: ± iij 


/o(1370)7r+ 


32.4 ± 7.7 ± 1.9 


30.4 ± 6.9 


34.7 ± 7.2 




0.76 ± 0.11 ± 0.03 


0.72 ± 0.11 


0.80 ± 0.11 




(198 ± 19 ± 27)° 


(156 ± 19)° 


(159 ± 14)° 


p°(1450)7r+ 


4.4 ± 2.1 ± 0.2 


5.8 ± 2.2 







0.28 ± 0.07 ± 0.01 


0.32 ± 0.06 







(162 ± 26 ± 17)° 


(144 ± 20)° 





"^/o{980) 


977 ± 3 ± 2 


976 ± 3 


974 ± 3 


(^ivie V / c J 








S'tt 


n no _L n ni i n m 
u.uy ± U.Ui ± U.Ui 


U.iU ± U.Ui 


n no 1 n m 
U.uy ± U.Ui 




f) no _|_ f) (14 _|_ f) no 


_n no _|_ f) OA 


OA 4- OA 


"^/o{1370) 


1434 ± 18 ± 9 


1401 ± 19 


1406 ± 15 


(MeV/c2) 








r/o{1370) 


172 ± 32 ± 6 


180 ± 34 


176 ± 30 


(MeV/c^) 










71.8/68 


93.0/68 


103.5/68 


C.L. 


35% 


2% 


0.4% 


-2 InZ^^jjj. 


-3204 


-3184 


-3172 



Table 1: Results of the Dg Dalitz plot fits. Fit A, all resonances are allowed (sys- 
tematic error follows the statistical). Fit ^ and Fit C do not include the p°(770)7r"'" 
and p°(1450)7r+ amplitudes, respectively. 
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Figure 3: S12 and S13 projections for data (dots) and Fit A (solid). The hashed 
area is the background distribution. 



ly-annihilation amphtudes or strong rescattering in the final state. In fact this decay 
is not observed in the —>■ K^K'tc^ final state[jlO|, pointing to the /o(1370) being 
a non-ss particle, as suggested by the naive quark model |T|]. 

In summary, Fermilab experiment E791 has measured the branching ratio of 
the decay Df 7r-7r+7r+ relative to Df 07r+ to be 0.245 ± 0.028tE];|jl^. We 
measure the mass and width of the /o(980) and /o(1370). Our results for the /o(980) 
parameters are = 0.09 ± 0.01 ± 0.01, qk =0.02 ± 0.04 ± 0.03, and mo = 977±3±2 
MeV/c^ Using the same Breit-Wigner function for the /o(980) as for the other 
resonances, we find Fq = 44 ± 2 ± 2 MeV/c^ and mo = 975 ± 3 MeV/c^. For 
the /o(1370) we find mo =1434 ± 18 ± 9 MeV/c^ and Fo = 173 ± 32 ± 6 MeV/c^. 
Finally, the fit of the Dalitz plot shows a dominant contribution from the /o(980)7r+, 
significant contributions from the /o(1370)7r"'" and /2(1270)7r"^, small contribution 
from the p°7r+ channels, and a negligible non-resonant component. The isoscalar 
plus TT"'" components correspond to over 90% of the 7r^7r+7r+ decay width. 
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